Quantum Multibaker Maps: Extreme Quantum Regime 



Daniel K. WojcikQ 
Institute for Physical Science and Technology, 
University of Maryland, 
College Park, MD, 20742, USA and 
Centrum Fizyki Teoretycznej Polskiej Akademii Nauk, 
Al. Lotmkow 32/46, 02-668 Warszawa, Poland; 



J. R. Dorfman 
Institute for Physical Science and Technology, 
and Department of Physics 
University of Maryland, College Park, MD, 20742, USA 
(Dated: February 1, 2008) 

We introduce a family of models for quantum mechanical, one-dimensional random walks, called 
quantum multibaker maps (QMB). These are Weyl quantizations of the classical multibaker models 
previously considered by Gaspard, Tasaki and others. Depending on the properties of the phases 
parametrizing the quantization, we consider only two classes of the QMB maps: uniform and ran- 
dom. Uniform QMB maps are characterized by phases which are the same in every unit cell of 
the multibaker chain. Random QMB maps have phases that vary randomly from unit cell to unit 
cell. The eigenstates in the former case are extended while in the latter they are localized. In the 
uniform case and for large h, analytic solutions can be obtained for the time dependent quantum 
states for periodic chains and for open chains with absorbing boundary conditions. Steady state 
solutions and the properties of the relaxation to a steady state for a uniform QMB chain in contact 
with "particle" reservoirs can also be described analytically. The analytical results are consistent 
with, and confirmed by, results obtained from numerical methods. We report here results for the 
deep quantum regime (large h) of the uniform QMB, as well as some results for the random QMB. 
We leave the moderate and small h results as well as further consideration of the other versions of 
the QMB for further publications. 

PACS numbers: 05.60.Gg, 05.60.-k, 03.65.-w 



I. INTRODUCTION 



The quantum mechanics of classically chaotic systems, 
often called quantum chaos, is by now a highly developed 
subject with an enormous literature, including mono- 
graphs by Gutzwiller [Q, Stockmann and Haake 
among others. The subject has been greatly advanced, 
as is usual, by detailed analyses of simple model sys- 
tems such as kicked rotors, quantum flows on surfaces 
of constant negative curvature. Harper models, and so 
on. Some of the central problems that have been stud- 
ied using these models include those of: (1) finding ex- 
planations for the efficacy of random matrix theories, 
(2) understanding the differences between quantum and 
classical transport, especially when Anderson localiza- 
tion plays a role in the quantum system, (3) studying 
the properties of quantum systems in the semi-classical 
limit, and (4) determining the role of decoherence in pro- 
ducing classically chaotic behavior of a quantum system 
as Planck's constant tends to zero. 

The present paper treats the quantum versions of sim- 
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pie model systems, multibaker maps, that have been used 
to study transport phenomena in classically chaotic sys- 
tems. A multibaker map consists of a chain of two- 
dimensional baker maps which are interconnected by 
means of a simple change in the baker dynamics. In the 
usual baker map on a unit square or torus, two vertical 
strips are stretched (by a factor of two) in the horizontal 
direction, contracted (by a factor of two) in the vertical 
direction, and the resulting horizontal strips are placed 
one above the other, in order to reconstruct the unit 
square. In the multibaker chain, each of the two hori- 
zontal rectangles are sent to adjacent cells, one to the 
right and the other to the left (Fig. |l|). Modified multi- 
baker chains have also been studied where there may be 
more strips and/or a more complicated dynamics includ- 
ing both area preserving and area non-preserving dynam- 
ics. These classical models provide simple, deterministic 
models of one dimensional random-walk processes with 
both diffusive transport and chaotic dynamics. They 
have been used to study connections between transport 
properties such as transport coefficients, and irreversible 
entropy production, and the chaotic properties of the 
models |, |, I 0, |, |, 0, |ll|, |l|, H [| @. 

Quantum versions of baker maps are well known and 
studied in some detail for a range of values of Planck's 
constant. Here we add a mechanism for transport of 
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FIG. 1: Classical multibaker model. 



probability amplitudes along a one dimensional chain of 
quantum baker maps. This quantum version of the multi- 
baker map (QMB) provides one realization of a quan- 
tum random walk process. In this paper we will con- 
centrate on the "most quantum" version of the QMB, 
obtained by using the largest possible value of Planck's 
constant h — l/2^n the Weyl quantization of the ordi- 
nary baker map . Our goal here will be to explore 
the transport properties of the QMB for this value of h, 
and for two different versions of the model, obtained by 
taking advantage of some phase-related arbitrariness in 
the quantization of the map. Later papers will explore 
further properties of QMB's including the semi-classical 
case, not considered here pi] . 

Area-preserving maps on torus admit a two-parameter 
family of Weyl quantizations , where the two 

parameters can be chosen to be phases. One can think 
of the two phases as offsets of lattice points that define 
the spatial and momentum coordinates of the map. If 
we choose the same phases in each unit cell of the chain, 
we have a "uniform" QMB. If we choose random phases 
from cell to cell, we obtain a "random" QMB. The differ- 
ent versions have quite different properties, as one might 
expect. The uniform QMB has many features in com- 
mon with those of continuous one-dimensional systems 
with periodic potential, including extended eigenstates, 
and ballistic transport, while the random case exhibits 
the usual phenomena associated with localization. Nev- 
ertheless, there are some interesting surprises, as we shall 
see in further sections, associated with transport in open 
systems. 

There are a number of formulations of quantum ran- 
dom walks already in the literature. We mention, in par- 
ticular, work of Aharonov et al. [Q, work of Godoy et 
al. and work of Barra and Gaspard |2^. The 

papers of Godoy and co-authors as well as that of Barra 
and Gaspard have interesting parallels with ours. These 
authors consider the motion of a quantum particle along 
a one-dimensional, periodic chain of scattering sites. The 
scattering sites are characterized by transmission and 
reflection amplitudes, which for a periodic system, are 
taken to be the same for each site. Godoy and co-workers 
consider the wave functions for their systems at discrete 
positions and discrete times, and propose a set of equa- 
tions similar to the ones considered here. These equa- 
tions are then solved using stationary phase approxima- 
tions, and the connections with Landauer's formula are 



discussed, for various parameters and particle statistics. 
Barra and Gaspard also consider a model similar to ours, 
and they analyze the scattering resonances for a finite, 
open system. By applying transfer and S matrices, they 
obtain expressions for the widths of resonances and the 
Wigner time delay, as functions of the system size. Their 
equations are in fact quite similar to ours, and a num- 
ber of results differ in the two cases only because of the 
differences in the details of the model studied. Their 
model has two channels per cell (particles moving to the 
left or right) but the particles can have a wide range of 
energies, and in some instances the high-energy limit is 
considered. The similarities with the work of Godoy et 
al. occur because the version of the quantum multibaker 
model considered in the present paper is the simplest 
possible, while more complicated versions, to be consid- 
ered in further papers have no direct counterparts in their 
work. 



Despite the similarities between our work and that of 
other authors, the focus of the work mentioned above 
generally differs from ours. We are particularly inter- 
ested in comparing and contrasting quantum and clas- 
sical multibaker maps, and in generalizing the QMB in 
a number of directions. These include an examination 
of the behavior of the QMB for smaller values of h in- 
cluding the semi-classical limit, and looking for traces, 
if any, of the chaotic classical behavior in the quantum 
version. The present paper is designed to identify impor- 
tant quantum phenomena that differ from those of the 
classical multibaker at large h, but which are expected 
to approach the classical results as the Planck's constant, 
h — 1/N, tends to zero. 

The plan of the paper is as follows: In Section |^ we 
will define the classical version of the model and study 
the evolution of piecewise constant probability densities. 
In Section III we quantize the multibaker map using Weyl 
quantization. There we will define the uniform and ran- 
dom QMB's and obtain expressions for the time depen- 
dent propagator appearing in the discrete time version 
of Schrodinger's equation. In Section IV we describe the 
behavior of the uniform QMB for h = 1/2. We find the 
eigenstates for both closed and open systems, as well as 
the steady state solutions for systems with particle reser- 
voirs at their boundaries. We then consider the transport 
properties of particles in these chains. Of particular inter- 
est in this connection is our finding that for open chains of 
uniform multibakers, and with absorbing boundary con- 
ditions, the escape of particles from the chain is sub- 
diffusive despite the ballistic transport of particles from 
the interior of the chain to its boundaries. We then turn 
to a brief discussion of the properties of random multi- 
bakers and show that the assumption of random phases 
leads to localized wavefunctions with very different prop- 
erties from the uniform case. Our results are summarized 



and discussed in Section VI 
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II. THE CLASSICAL MULTIBAKER MAP 



A. Closed, periodic case 



The classical multibaker map provides a reversible, de- 
terministic realization of a one-dimensional random walk. 
It is the simplest area-preserving, deterministic model 
for diffusion of a particle on a one dimensional lattice, 
whereby the particle makes steps either to the right or 
left at equally spaced time intervals. The multibaker 
map can be adjusted for any set of step probabilities, 
p,q — 1 — p,0 < p,q < 1, where p is the probability 
of making a step to the right. The classical multibaker 
map is based upon the usual baker's map, B on the unit 
square, {0 < x,y < 1), defined by 



B(x,y) 



{x/p,py), 

{{x-p)/q,p + qy), 



for < X < p, 
for p < X < 1. 



(1) 



The multibaker map is constructed by taking a linear 
chain of L adjacent unit squares, labeled by the index 
n, n = 0, . . . , L — 1, such that any point on the chain is 
labeled by the three quantities, n, x, y with < x,y < 1. 
Then the action of the map, M on any point is obtained 
by combining a baker's map with translation of each rect- 
angle to the right or left, as given by 

M{n,x,y) — {n + \,x/p,py), ior < X < p, (2) 
= {n — \,{x ~ p) I q,p + qy), ioY p < X < 1. 

This arrangement has the property that there is a prob- 
ability p of choosing a point which moves one square to 
the right, and probability q of choosing a point which 
moves one square to the left. To complete the speci- 
fication of the map, one must append boundary condi- 
tions to the transformation given by Eq. (||). Such condi- 
tions may include periodic, or absorbing boundary con- 
ditions, or one might specify that the ends of the chain 
are connected to reservoirs which maintain a constant 
density of points at the boundaries. As a chaotic sys- 
tem, the multibaker map is a measure preserving map 
with positive and negative Lyapunov exponents, given 
by \± = ±[pln(l/p) + gln(l/(7)]. This map has been 
used to study the properties of deterministic diffusion in 
a chaotic system, studies of the connection between dif- 
fusion coefficients and Lyapunov exponents for an open 
chain, a study of entropy production in the relaxation 
to a uniform equilibrium state, and has been extended 
to provide simple models for viscous and heat flows as 
weU |, |, I 0, 11 0, |ll| [1| Q @. 

The classical version of the quantum multibaker con- 
sidered here was discussed in |^ . We consider here a sim- 
ple form of this classical model constructed to be a classi- 
cal version of the h — 1/2 quantum system. We will study 
the evolution of probability densities integrated along the 
stable direction (y) and piecewise constant on two halves 
of every multibaker cell (along the unstable direction). 
This space of densities is 2L dimensional. Therefore, the 
evolution operator for this class of probability densities 
has the same dimension as the quantum multibaker prop- 



We consider the classical evolution of phase space den- 
sities under the dynamics given by Eq. (0) with p = 1/2. 
Since the quantum version will describe probability am- 
plitudes in either space or momentum, the classical coun- 
terparts are obtained by projecting the classical densities 
along the x or y-axes, respectively. We restrict our atten- 
tion to probability densities projected onto the unstable 
x-direction and we take them to be constant on intervals 
< a; < 1/2, 1/2 < a; < 1, n = const, to mimic the 
h = 1/2 quantum case. Then the projected distribution 
is 



g{n,x,t) := / g{n,x,y,t)dy 



gL{n,t), 
gR{n,t), 



for < a; < 1/2, 
for 1/2 < x < 1, 



(3) 



and it satisfies a Frobenius-Perron equation given by 
gL,R{n,t + ^) = \[gL{n-l,t) + £»i^(r^ + l,t)], (4) 

with periodic boundary conditions gL.Ri'n- + L,t) = 
gL,f({n,t). Since this equation is linear we may suppose 
that gL.Rii^it) represent the deviations from a uniform 
equilibrium state, and may take both positive and nega- 
tive values. An eigenstate of the right hand side of Eq. 
(^) satisfies 

>^gL{n) = XgR{n) = ^[gL{n - 1) + + 1)]- (5) 

It follows that either A = or gL{n) ~ gR{n). Clearly, 
the L vectors of the form 

gdk-l) = -gR{k + l)^Q, 
gL{n^k-l) = gB{n^k + l) = Q, 

belong to the kernel, A = 0, of the classical discrete 
multibaker. For the case where A 7^ we look for so- 
lutions of the form gL,R{n) = Ae^"^^ and Eq. leads 
to A = cos??. The general solution is, then, gL,R{n) = 
Ai cos(i?n) + A2 sin(i9n), where periodic boundary condi- 
tions lead to d — 2kTr/L, and the normalized eigenstates 
can readily be determined. For odd L = 2M -|- 1 we have 
the L solutions 



1. M solutions of the form gL,R{n) — A cos 2^ 
cos2|iL; k^l,...,^] 

2. M solutions of the form gL,R{n) = Asm 
cos^; l,...,i^; 

3. 1 solution gL.R.in) = A, A = 1; fc = 0. 
For even L = 2AI we have the L solutions 

1. M — 1 solutions of the form gL.Rin) = A cos 



A = 



2k7:n 



agator considered in Section III 



A 



L ' 



L 
■ ' 2 



1; 
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2. M - 1 solutions of the form QL,R{n) = Asin^^, 
A = cos^; fc = 1,. - 1; 

3. 1 solution QL,R{n) = ^, A = 1; fc = 0, and, 

4. 1 solution QLM{n) = (-1)"^, A = -1; fc = Af = 
L/2. 

We see that in the odd case there is an approach to equi- 
librium: all the eigenvalues have absolute value strictly 
less than 1, apart from the one corresponding to the 
uniform distribution. The even case is sensitive to the 
"even-odd" oscillations of the location of a point along 



the chain. These oscillations can be removed by combin- 
ing two successive steps. 

B. Open case (absorbing boundary conditions) 

For the open chain with absorbing boundary condi- 
tions, the dynamics inside is the same as in the closed 
case and is given by Eq. (H), therefore the general solu- 
tion is also given by QL,R{n) = Ai cos{'dn) + A2 s\n{dn). 
Absorbing boundary conditions Qr^l[—^) — Qr,l{L) = 
lead to 



X[Ai cos(L - + A2 sm{L - 1) 



- [Ai cos 1^-1-^2 sin 1?] , 

^[Ai cos(i - 2)?? -I- A2 sin(L - 2)i9], 



(6) 
(7) 



where A — cos-d. They have nontrivial solutions if and 
only if sin(iy + Ij-d = 0, leading to i?fc = -^fj, where k — 
-L,...,-1,1,...,L, and A(-79) = A(i?). Thus finally 
k = 1, . . . ,L, which gives the L states of the form 



QL,R{n) = A sin 



k{n + l)7r 
L + 1 



(8) 



The remaining L states are in the kernel: L — 2 are of the 
form given by Eq. with n = 1,...,L — 2. Two other 
states corresponding to A = are ql{L — 1) = 1, PL(fc 7^ 
L-1) = Quik) = 0, and ^^(0) = 1, QR{k ^ 0) = QL{k) = 
0. Thus we easily obtain a spectral decomposition for the 
simple operator treated here, with absorbing boundary 
conditions. 

The probability of finding the particle in the system 
decays with the escape rate 
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lim 

t — >oo 



log Pit) 



t 



(9) 



where P{t) :— g{n, t), given by the largest eigenvalue 

^2 



7 ^ 



log|cos— — I » — 



1' 



2Z2' 



(10) 



for large L. 



The open, discrete multibaker with reservoirs 



where gi , g2 are the incoming densities of the left and 
right reservoirs, respectively. A solution is found imme- 
diately by observing that in the steady state ghin) = 
gR{n) = gin) and that Eq. (^ leads to gin +1) — 
2gin) — gin — 1). A solution satisfying the boundary con- 
ditions is therefore 



gin) 



g2 + Lgi nig2 - gi) 



(13) 



This linear profile expected from the Pick's law ^, ||] 
is shown in Figure B. 
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FIG. 2: Linear profile of the classical discrete multibaker map 
of length L = 101 with reservoirs (qi — 0.1, f>2 = 0.9). The 
horizontal axis range is [—5, 105]. 



Next we connect particle reservoirs to a finite chain and 
look for steady state solutions. These are time invariant 
solutions to Eq. (H), with the boundary conditions 



eUO) = gRiO) ^ ligi + gRii)h 



(11) 



QLiL-l) 



gRiL-l)^-[gLiL-l) + g2], (12) 



III. QUANTUM MULTIBAKER MAP: THE 
GENERAL MODEL 



In order to quantize the multibaker map, we start with 
the quantum baker map, as described by Balazs, Voros 
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|l9[ , and Saraceno IQ, and then produce a quantum 
multibaker map by forming a chain of unit squares, ap- 
plying the quantum baker map to each square, but trans- 
ferring the new quantum states to the adjacent squares 
according to the procedure used in the classical case. 

The method for constructing a quantum version of 
the regular baker map is as follows. We consider the 
x-direction to be the "spatial" direction of the system, 
and the y-direction to be the "momentum" direction. 
Then the number of quantum states N in the unit square 
should satisfy N — PQ/2ttH, where Q = 1 is the spatial 
extent of the unit square, and P = 1 is the range of mo- 
menta. This leads to the simple formula, h = 1/{27tN), 
where N is an integer. We usually, but not always, 
take TV to be an even integer, so that one half of the 
quantum states can be associated with each half of the 
unit square. One then constructs a set of N "posi- 
tion" states for a unit square, with position eigenvalues 
q, = 2T:h{j + ifq) = {j + ^q)/N, j = 0, . . . , iV - 1 and 
a set of N "momentum" states with momentum eigen- 
value pk = 2iTh{k + ipp) ^ {k + Lpp)/N, k = 0, . . . ,N - 1. 
We require that < ipq,p < 1. The position and momen- 
tum states are related to each other by means of a simple 
Fourier representation, with N terms, given by 



space 



{GN)kj =< Pk\qj >= 



1 



-27ri(fc+Vp)(j+V,)/A^ 



N 



(14) 



We include subscripts on p, q in the notation for the 
Dirac matrix element to identify the integers which are 
attached to the p and q representations. The phases, 
fq,fp are as yet unspecified. In the literature on the 
quantum baker map, these phases are often taken to be 
(the simplest [|l9|) or 1/2 (most symmetric map [pO[). 
Here we will take advantage of the possibility to choose 
these phases so as to represent different situations that 
may have some relevance to physical phenomena. 

The time dependence of the quantum baker map is de- 
termined by constructing a propagator for the change in 
the quantum states over one time step. This propaga- 
tor consists of two parts: First one transforms the "left" 
part of the Hilbert space (in the position representation) 
into "bottom" subspace (in momentum representation) 
and the "right (position) part" into the "top (momen- 
tum) part" . Then one uses the Fourier relation between 
position and momentum states, Eq. (^^, to change the 
basis from momentum back to position representation. 
The first transformation consists of two Fourier trans- 
forms over A^/2-dimensional space, the other is the in- 
verse Fourier transform over the whole A^-dimensional 



Gn/2 

Gn/2 



(15) 



The action of B on a position-space wave function is 
understood as follows. We represent the position-space 
function as a column vector with N elements, the top 
N/2 elements referring to the quantum states with num- 
bers j = 0, 1, . . . , {N/2) — 1, which we denote as the "left" 
states. The bottom N/2 elements having quantum num- 
bers j = N/2, N/2 1, . . . , iV - 1, are called the "right" 
states. The block diagonal matrix, with blocks Gn/2 ap- 
pearing on the right hand side of Eq. (^5|) transfers the 
left and right spatial states to the "bottom" and "top" 
momentum states, respectively, according to 



>6(t+l)' 









GN/2{Vq,Vp) 

GN/2iVq,Vp) 



*r(t) 

(16) 

This operation defines the quantum baker map. How- 
ever, we are left with a quantum state in the momentum 
representation. We now change the momentum state rep- 
resentation into a spatial state by means of the matrix 
asinEq. (0). 
Finally we can construct a quantum multibaker map 
(QMB) by considering a chain of unit squares, each taken 
to be an individual quantum system, but which exchange 
quantum states according to the rules of the quantum 
baker with an interlacing process formed in analogy with 
the classical multibaker map Eq. (||). That is, the po- 
sition space functions at site n are transformed to mo- 
mentum space functions at sites nil, according to the 
rule 



fb{n+l,t+l) 
^t{n-l,t+l) 



GN/2{n) 
GN/2{n) 



Here we can allow for the phases (pq^^ 



*/(n,i) 

*r(«, t) 

(17) 

to vary from 



one cell, denoted by n, to the next, and we incorpo- 
rate them in the transformation operators GAr(n) = 
G'Ar(<y9g(n), (^p(n)) at that site. After this transformation 
is carried out, we change from the momentum to the po- 
sition representation at each site according to the same 
rule as in an ordinary quantum baker map, that is 



■^i(n,t- 
^r{n,t 



= G-^\n)- 



fb{n,t+l) 
■^t{n,t+l) 



(18) 



J 



Thus in the position representation the quantum multibaker map is given by 



■^i{n,t + l) 
*r(n,i + 1) 



= GN\n) 



GN/2{n-l) 

Gjv/2(" + l) 



Explicitly, we have 



*,(n - l,t) 
■^r{n + l,t) 



i + 1) = [GN\n)]i^b ■ GN/2{n - 1) • *i(n - 1, t) + [G^^n)];,* • Gw/2(n -I- 1) • ^rin + 1, t). 



(19) 
(20) 
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t + 1) = [Gj,\n)]r,b ■ GN/2{n -l)-^i{n- 1, t) + [Gj,\n)]r,t ■ Gjv/2(" + 1) • + 1, t). 

I 



(21) 



Here, in an obvious notation, the matrices [G^^Jq^/j are 
N/2 X N/2 block sub- matrices that comprise G^^. The 
general case can be treated numerically, of course, once 
the phases are specified. It is of interest to consider the 
special case N = 2, since much of the work can be done 
using simple analytical methods, and since this case cor- 
responds to the largest possible value for Planck's con- 
stant, namely, h = 1/2. This is the case we study here. 

The local dynamics are characterized by the two 
phases, i^q,ipp, which parameterize the Weyl quantiza- 
tions of the baker map. If we take the same pair of phases 
at each site we obtain the uniform model. If we choose 
them randomly from some distribution at each of the 
sites, we get the random model. In this paper, when we 
treat the random model we will assume that the phases 
are chosen according to a uniform distribution on the 
unit circle Q. 

A complete specification of the model is obtained by 
adding the boundary conditions to the above equations. 
In this work we restrict our attention to the closed case 
(with periodic boundary conditions), and open cases 
(with either absorbing boundary conditions or with "par- 
ticle" reservoirs at the ends of the chain). 



IV. UNIFORM QUANTUM MULTIBAKER 

The uniform quantum multibaker is characterized by 
a set of phases ipq , ipp that are independent of the site 
index, that is, they are the same for each of the trans- 
formation matrices generating the map, as described in 
Eq. (p^). This makes transport in the uniform multi- 
baker chain similar in many respects to transport in a 
one-dimensional periodic solid. Here we solve this model 
for time dependent and stationary quantum states with 
appropriate boundary conditions: the closed, periodic 
chain; the open chain with absorbing boundary condi- 
tions; and the open chain attached to leads at each end, 
producing a stationary, non-equilibrium state. We begin 
with the periodic chain. 



A. Closed, periodic case 

We consider the periodic, uniform multibaker chain, 
with L sites and N = 2. The equation connecting the 
quantum states at time t + 1 to those at time t is 

*z(n,t + l) = /o(.goo*i(n- 1,0 +5oi*r(ri + l,t)), 
^rin,t + l) = /o(5io*i(n-l,t)+5ii*r(ri + l,t))- 



with (0 < (pq,(pp < 1), and fo,gki given by 
/o = (Gi(^„^p))oo = 

1 



gkl = (G2 {ipq,ipp))kl 



T2 



(22) 



Since the system is periodic, Bloch's theorem guarantees 
the existence of eigenstates of the form 



(23) 



Periodic boundary conditions, '^{L) = \E'(0), imply that 
■d — 2k'iT/L^ k = 0, . . . , L — 1. Clearly, A is an eigenvalue 
of the quantum multibaker propagator if and only if 



500 /oe - A .901 /oe"' 
510/oe-''' Sii/oe^" ~ A 



Using the notation 

a = {l+ipq + (y9p)7r/2, 

P = {I + iPq + (Pp - 2(pqipp)TT/2 - 

we find that 



0. 



A 



T2' 



[cos(i9 + a)±iJl + sin^(i? + a)]. 



(24) 



(25) 



(26) 



Note that (24) can also be written as 



A/e 



1/3 



7A 



v + - = 

V 



1 

V2 
1 



1 

U+- 

u 



(27) 



where u ~ xe", v ~ A/e*^. Since ■& is real it follows that 
|A|2 = 1, so A = e^^ and 7-/3 £ [7r/4, 37r/4]U[57r/4, 77r/4]. 
Therefore the "quasi-energies" , 7, lie in two bands of 
length 7r/2 symmetric with respect to the center of the 
unit circle. The exact location depends upon the phases 
ipq,ipp. Making use of the boundary conditions we obtain 
the eigenvalues of the closed multibaker map 

^ T2 



cos(a + 2kTT/L) ±1^/1 + sin^(a + 2kTT / L) 



(28) 

with a and (3 given by Eq. (|25|). The correspond- 
ing eigenstates are given by Eq. (^3|) with the con- 
stants connected by Ar — A;[sin(a -I- 2fc7r/L) =F 

^ 1 -I- s\n^{a + 2fc7r/L)]e'('^('^''-'^«)/2-2fc'^/-E') . 

When a is an integer multiple of tt/L the spectrum 
is doubly degenerate. This non-generic case happens for 
instance for the most common choices of phases ((/3g — 
(y9p = or 1/2). The quasi-energy spectrum of the closed 
uniform multibaker for the phases Lpq = Lpp = 1/2 is 
shown in Figure |^. 
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FIG. 3: Eigenspectrum of the closed uniform quantum multi- 
baker for the chain of length L = 101 cells with periodic 
boundary conditions. Phases are ifq=tfp = 1/2, a = n, [3 = 
2,-kIA. 



B. Open case: absorbing boundary conditions 

Next we consider the uniform quantum multibaker, 
still for N = 2, but with open boundaries. In the classi- 
cal case, open boundaries are important for the applica- 
tion of the escape-rate formalism of Gaspard and Nico- 
lis which relates the rate of decay of the initial num- 
ber of particles on a large, open chain to the diffusion 
coefficient, and then to the Lyapunov exponents and the 
Kolmogorov-Sinai entropy of trajectories on a fractal re- 
peller, i.e. the set of initial points for trajectories that 
never leave the chain It is of some interest, then, 

to contrast the classical and quantum cases. 

We take the multibaker dynamics given by in the 
cells n = 1, 2, . . . , L — 2. At the boundary cells we allow 
the probability density to escape from the right half cell 
for n = 0, and from the left half cell for n = L — 1, and 
nothing enters the system from the outside. The latter 
condition requires 



*b(0,t) = 0, 
^t{L-l,t) = 0. 



(29) 
(30) 



Due to the escape of probability density, the eigenval- 
ues that determine the time dependence of the probabil- 
ity density in each cell move to the interior of the unit 
circle. A simple proof of this fact is given in Appendix 
We show there also that the kernel is two-dimensional. 

To determine the non-zero eigenvalues < |A| < 1 of 
the open chain, we first write the eigenvalue equation in 
the momentum representation. Then every eigenstate ^ 
satisfies the equation 



" *b(n + 1) ■ 




7o/A 





. *t("-l). 







/o/A 



■G2 \Vq,Vp) 



(31) 

Viewed in terms of the "top" and "bottom" states, we 
see that the solution of Eq. (|l]) can be neatly formu- 
lated as a scattering problem (see Figure ^), where the 
incoming waves are "^bin) and "iltin), and the outgoing 




incommg outgoing 
FIG. 4: Scattering from one cell. 



waves are 5'(,(n + 1) and ^'t(n — 1) with a one-cell scat- 
tering S'-matrix, and a one-cell transfer T-matrix. Those 
are defined, respectively, by the relations 



" *4(n- 


!)■ 


= 5- 


■ *f,(n) " 




1)J 





and 



"if bin + 1) 



*t(n- 1) 



(32) 



(33) 



Explicitly, the S'-matrix is, for the uniform multibaker, 
given 



72^ 



/A 



and the T-matrix is 



T = 



(34) 



(35) 



We find it convenient to use the transfer operators, T, 
to carry out the determination of the eigenvalues. A, gov- 
erning the rate of decay for an open system. To do this 
we first use the transfer operators to relate the quantum 
states at one end of the chain to the states at the other 
end, and then use the open, absorbing boundary con- 
ditions to obtain an explicit equation for A. First, the 
states at the two ends of the chain are related by 



fb{L-l) 
*t(i-2) 



T 



L-2 



*t(0) 



(36) 



To use the boundary conditions, we first look at cell L — 1. 
We note that in the open multibaker (L — 1) = 0. Then 
using Eq. (|l]), we obtain $f(L- 2) = (/o/A) [gio$b(i- 
1) -f gii$t(L - 1)], and ^t{L - 1) = [(A//o)*t(L - 2) - 
gw^b{L - l)]/gii. Thus 



= [ -.910 (A//o) 



(37) 



Using Eq. ( |32| ) and (|3q), we easily find that ^'b(l) = 
(/o5oi/A)^t(0). Thus the equation that determines the 
decay rates is 



= [ -/o5io A ] • T 



L-2 



A 



(38) 
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where a scalar product of matrices is to be taken as indi- 
cated. To get a useful form for this equation we need to 
find the eigenvalues of the transfer matrix T. We denote 
the eigenvalues of T by x+tX-j which are obtained as 
solutions of the quadratic equation 



1 



(39) 



where a,P are given by Eq. (p5|). Using, as before, the 
notation u = xe*"; v = A/e'^ we obtain the same formal 
relation between u and v as in the periodic case, Eq. (27), 
i.e. 



1 

71 



(40) 



The two solutions m+,u_ satisfy u+u- = 1, and m+ + 
M_ = %/2[w + l/v]. Since \v\ = \X\ < 1, it follows that 
do not lie on the unit circle. In particular, they 
must be different and so the matrix T is non-degenerate. 
We take \u+\ > 1 > to define them uniquely, and 
use u± = x±e*". If we set u± — e****, and then solve for 
V we obtain 



v± — —= [cos K ± I \/ 1 + sin^ n\ . 
V2 



(41) 



Interesting solutions are those where k is not purely real, 
that is, K S C \ R. We next use a simple identity for the 
L-th power of non-degenerate matrix T, given by 



xl 



X- 



x+ - X- 

g-ia(L-l) 



■T 



x-x\ 



X+X- 



to write Eq. 



2v' 



2 + ^ 



sm K 
is]) in the form 

sin(L — 2)k = 



X+ - X- 
[sin(LK)T-sin((L- 1)k)/], 



V2tJ 



\/2t; 



where v is one of v±. With the help of 
can reduce (E2|) to 



sin(L — 3)k, 
(42) 

and (41) we 



siniK-|-sinK cos(L— V)K^ie\/ \ + sin^ Ksin(i — 1)k = 0, 

(43) 

which can be further reduced to 



,2L 



— 1 — 2sinK[sinK — eVT+sii?K] = 0. 



(44) 



In the above equations e = ±1 corresponds to the sign 
in (^). Clearly, we can get all of the possible solutions 
multiplying Eq. (44) for two different signs. This leads 
to a very simple equation 



sin^ Lk 



sin K = 0. 



(45) 



The only real solutions of this equation are k = fevr, k G 
Z, but, as mentioned above, they must be discarded. If 
we write Eq. (|45|) as 



sin Lk = 1(5 sin k, 



(46) 



where b = ±1, we can treat it as a "perturbation" in 
5 of equation sini^ — [Q. Thus we can obtain the 
solutions of interest by expanding k in powers of 5 about 
the values k = k-n/L, k = 1, ... ,L — 1, |^5| and then 
at the end, setting 5 = ±1. This approach gives results 
which quickly converge numerically, for all allowed values 
of fc. To apply this procedure it is convenient to rewrite 
Eq. mm in a polynomial representation. 



,2L 



1 - 8iu^(u - l/u) = 0. 



(47) 



Then, by taking u = exp(i(fc7r/L -I- 6ai + 5^a2 + . ■ .)), 
one can determine the coefficients a^, and check the con- 
vergence of the series numerically. Figure |^ shows the 
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FIG. 5: The first and the fifth order approximate solutions 
and the numerical solutions of Eq. ( |46| ) for L = 101. We took 
here 5 = 1. The numerical solution is the middle curve. Fifth 
order approximate solution is the one closer to the numerical 
solution. 

absolute values of the approximate solutions (in the first 
and fifth order) and the numerical solutions for L = 101 
sorted according to increasing amplitude. The first few 
coefficients in the expansion of u are 

ai = —b/L, 

a2 — —iab/L^, 

as = b^/6L + b(2-3b'^)/2L^, 

04 = ii2ab^/3L^ - ab{8b^ -3)/3L-^), 

where 



a = cos(A:7r/L), 



sin(A;7r/L). 



(48) 



Numerical studies show rapid convergence of amplitudes 
and slower convergence of phases. 

Next we calculate the approximate eigenvalues of the 
open quantum multibaker. Using A = e'^u, keeping < 
1 solutions, to second order in 5 we obtain 



9 



l^a + ie + b'^ ) exp 



I 



cxp ■ 



ieab"^ (2 &2 + 3) 
2i2 (1 + 62)^/2 



(49) 



where a, 6 are given by ([48|), while e = ±1 enumerates the 
solutions. The non-exponential factor on the right hand 
side is the unperturbed solution. Figure IV B shows the 
absolute value of v (in the fourth order approximation) . 
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the eigenvalue corresponding to fc = 1: 



V 



-3-2-10123 
argument (k) 

FIG. 6: The amplitude versus argument of A in the fourth 
order of approximation for the chain of length L = 101. 

The general result and the numerical studies suggest 
that the value of v that is the nearest to the unit circle 
occurs when k takes on one of the four values, k G {1, L — 
l,i+l,2i— 1}. To get the leading term in the large-size 
limit we find it convenient to use an alternative expansion 
of Eq. (^6|), in powers of L~^. While this leads to an 
asymptotic solution for v which quickly diverges for most 
fc, it gives us the correct leading order behavior for k <ti 
L. Then the expansion of the solution of Eq. ( ^6| ) in 
powers of 1/L, as Kk = fc7r/_L + + h^/ + . . . for 

small k yields 

Kk ~ kn{l/L + iS/L^ i5{l + k'^T:^/3)/L* + ...), 

(50) 

(51) 
(52) 



which gives the asymptotic formula for v{k) 



v{k) 



—= [cos Kk^i 
V2 



1 + sin^(Kfe)] 



exp 



±i 



2i2 



-fc2^2 



cxp ■ 



L3 



The escape of probability density from an open system 
asymptotically is dominated by the eigenvalue closest to 
the unit circle. Therefore the escape rate 



7 := — lim 

t — >oo 



log P{t) 



t 



(53) 



7 



-log|<l)p 



27r2 

L3 



(54) 



This result means that even though the motion in- 
side the quantum niultibaker is faster (ballistic) than in 
the corresponding classical system (diffusive), the ejfu- 
sion (decay of probability density) is slower than that 
for t he co rresponding classical system, ( Q| ; see also Sec- 
tion |nB| ) 



Tclass 



2i2 



(55) 



of the uniform quantum multibaker map is obtained from 



It is interesting to compare this result with those ob- 
tained by Barra and Gaspard |27| in their study of scat- 
tering resonances for an open, periodic chain of scatter- 
ers. In high energy limit they found that the logarithms 
of the magnitudes of the eigenvalues can be bounded 
above and below by functions that scale as 1/L. They 
expect that the lower bound, given by the eigenvalues in 
the middle of the band, should hold also for lower ener- 
gies Q. On the other hand, the upper bound, which 
gives the escape rate, is given by the resonances near 
the edges of the bands which are harder to estimate at 
low energies. Therefore this bound is more difficult to 
control. 

This reasoning is consistent with our findings. In our 
case, the eigenvalues of the smallest magnitude are those 
for whi ch to k sa ±L/2 (the middle of the band; see 
Figure IV J^ ). Thus their magnitude can be estimated 
from ([49|) setting a = 0,b = 1 and therefore their loga- 
rithms scale as 1/L. On the other hand, the eigenvalues 
of largest magnitude, which give the escape rate, lie at 
the edges of the band. 

The discrepancy between our results is not surprising 
for the high-energy limit corresponds to semi-classical 
limit for our system, and in the present work we con- 
sider the extreme quantum case. 



C. Steady state solution 

Suppose now that the multibaker of length L is con- 
nected at both ends to infinitely conducting leads. We 
suppose that there can be traveling waves in the leads 
moving to the right and to the left. These waves are 
most conveniently described in terms of the momentum 
space representation of the wave functions, and we re- 
call that the "bottom" states come from the left and the 
"top" states come from the right. Thus to the left of the 
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chain, n < 0, we take the traveling waves to be 

*b(n,t) = ^e*('^*-'="), (56) 
^t{n,t) = Se'M+M^ (57) 

and to the right of the chain, n > L — 1, the waves are 
^t{n,t) = c-e^M-fc")^ (58) 
^t{n,t) = De'M+fc"). (59) 

Here A, D are the ampUtudes of the incoming waves, 
while B,C are the amplitudes of the outgoing waves. 
Due to the dynamics on the multibaker, we can match 
the incoming wave functions in the leads to the proper 
momentum space functions for the unit cells at and at 
L — 1. This matching condition is simply 



where the elements of S'-matrix are 



*f,(0,t) = Ae'"* 
$t(L-l,t) = De 



(60) 
(61) 



We will use a scattering approach to find the outgoing 
amplitudes, _B, C for the steady state solution, as well 
as to solve the problem of the relaxation of some initial 
state to a steady state. First we consider the steady state 
solution for the baker chain with conducting leads. 

The steady state solution is defined by the condition 
that the time dependence of the wave function can be 
incorporated in a time dependent phase factor. Since 
the dynamics takes place at discrete times, there is a 
u) such that 'ift.birijt + 1) = e"^5't_h(n, <) which im- 



plies that '^t^i){n,t) = 
Eq.^, ^t{0,t) = 
ing *b(n,t) = 
equation 



4'f,fc(n, 0). In particular, using 
'^-6(0,0), so that 5 EE w. Writ- 
*'i'b{n-), etc., we obtain the steady state 



" + 1) ■ 









foe 






Joe 



(62) 



The transmission and reflection coefficients for the 
chain can be expressed in terms of the scattering S'- 
matrix, given by 



%(-!) 
*b(i) 



*6(0) 

L*t(i-i) 



(63) 



So,L-l 



ro,L-i t'o,L-l 



' O.L-1 



(64) 



Here the unprimed coefficients refer to waves incident on 
the left end of the chain, while the primed quantities 
refer to the waves incident on the right side of the chain. 
The transmission and reflection coefficients, T,T' , R, R' 
respectively, are then obtained from the elements of S by 



T : 



|io.L-i| 



R^\ro,L-i\^ 



(65) 



and similarly for the primed quantities. Unitarity of S 
implies T = T' , R = R' . In order to calculate the S'- 
matrix, So,L-i, for the chain, we proceed as for the ab- 
sorbing case, by looking at the transfer and scattering 
matrices for one cell, and building up the matrices for 
the chain by iteration, cell by cell. Consider the cell la- 
belled by the index n. The S'-matrix for the n-th cell is 
given by 



Mn-l) 
*b(n + 1) 



Sri 



and the transfer T-matrix is 



*b(n + 1) 



^*b(n) 
*t(n- 1) 



(66) 



(67) 



Each of the matrices can be given in terms of the other, 
thus, 



, (68) 
, (69) 



The S and T matrices can easily be obtained by trans- 
forming the dynamical equations ( |6^ ) to momentum rep- 
resentation dlq), so that 



s = 


r 
t 


t' ' 
r' 


=> T = 


-t'-\ t'-' _ 


T = 


a 

J 


7 
S 


^ S = 


-5-^(3 ' 



*t(n-l) 



foe~ 



900 goi 
9w 911 



Jb(n) 
*t(n) 



(70) 
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from which the ^-matrix follows as 



V2 



jTripg{l-ipp) 



(71) 



The T-matrix is then given by 



(72) 

The scattering matrix for the whole multibaker 5o.l-i 
can easily be derived from Tq^l-i := ■ ■ ■ ■ ■ Ti ■ Tq. 

Its unitarity can also be verified. For the uniform system 



X4 



X- 



■T 



X-Xa 



X+X- 



x+ - X- x+ - X- 

where x± ^re roots of characteristic polynomial of T 



(73) 



X± 



= e-'" [ V2 cos(/3 - w) ± Vcos2(/3-cj)] , (74) 



and a,/3 are given by Eq. (|25|). Depending on the sign 
of cos2(/3 + w) there are two types of solutions: if the 
frequency of the incident wave falls in one of the quasi- 
energy bands 

cos2(/3-cj) <0^cj-/3e [7r/4,37r/4] U [57r/4,77r/4], 

(75) 

we have the oscillatory case with some interesting struc- 
ture. Otherwise, when the frequency of the incident wave 
falls in the gap, we observe almost total reflection of par- 
ticles coming from the leads to the chain, becoming total 
as L cxD (the exponential case). 

1. If cos2(/3 — w) < (oscillatory case), the charac- 
teristic roots are: 



X± = e~'"[\/2cos(^ -Lo)± V-cos2(/3-tj)], (76) 



thus \x±? = 1. Set x± 
scattering matrix for the chain becomes 



So,L-l — 

ZL 



Sinn e 



— iaL 



sinK e*"^ 
siniKe'("~'''^j') 



To simplify the formulas we introduce 



(77) 



\/2sinnKe"'(''~") - sinK(n - 1) 
cos nn sin k — ie sin nn \/ 1 + sin^ n 



where £ = ± is the sign of sin(/? — oj). Then the 
transmission and reflection coefficients are 



R 



T 



sin^ Lk 



1 



sin^ K + sin^ Lk 1 + ' 
sm K 1 



sin^ K + sin^ Lk 1 



sin- Lk 



(78) 



(79) 



Some interesting special cases occur when: 

(a) K = kn + tt/2, L odd: T = 1/2; 

(b) K = fcTT + 7r/2, L even: T = 1 ; 

(c) K = kn: T = 1/(1 + ^2) ; 

(d) K = kn/L: T = 1. 

We will refer to the cases when T = 1 as transmis- 
sion resonances. They occur when sin Lk — 0. On 
the other hand, one can see from Eq. ( [TT] ) that the 
S-matrix has poles when Eq. ( ^5| ) is satisfied. Hence 
the poles of the S-matrix determine the eigenstates 
of open system. 

2. In the exponential case, when cos 2(/3 — w) > 0, we 
have |x±P > 1, andx+xl = 1, so that |x-| = j^- 
Then the transmission and reflection coefficients 



R 



T 



(Ix+l^-lx-l^)^ r«n^ 

(Ix+l^-lx-n^ + dx+l-lx-l)^ ^ ' 

l-|;^_|2(i-l)«l (81) 

(Ix+l-lx-l)^ r82^ 

(lx+l^-|x-l^)^ + (lx+|-|x-l)2 ^ ' 

|x-|'(^-^)«0. (83) 



D. Density profile in tlie steady state — violation 
of the Ficlc's law 



As mentioned above, the oscillatory case provides some 
interesting structures, illustrating the interference be- 
tween waves traveling to the right and left along the 
chain. The algebra is tedious but straightforward, and 
we don't reproduce it here, merely stating the final re- 
sults. 

The wave function in the steady state is 



*t(n) 



g— ia(n+l — L) 



[zL-n-i'biO) ~ e"^e'("-"'^'') sinnK ^t{L - 1)], 

sin K{n + \-L) l'b(O) -I- z„+i^'t(L - 1) 



rg-MLgi(a-7rVg) , 



ZL 
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We introduce the probability densities, ql and qr 
from the left and right leads, respectively, in terms of 
the corresponding wave functions, by writing \['b(0) = 

'^t{L — 1) = yg^e"*, where 77 denotes a relative 



phase between the wave functions at the two ends. Then, 
introducing the angle Lp = Tr{ipq — (pp)/2 + aL + r/, we ob- 
tain the total probability density at cell n 



, , sin^fL — n ~ 1)k + sin^fL — n)K, + sin^ k sin^ kit, + sin^ Kin + 1) + sin^ k 
Q{n) = ^ -^-^ QL + QR 



.y/QLQR 



{sin(L — 1 — n)K [zn+ie'^'^ — z^+iC *'^] — sinnK [z2_„e 



At resonance (n 



g{n) = 1 + 



: kiT / L) it takes form 
sin^ Kn + sin^ K(n + 1) 



sin^ K 



{ql + Qr) 



-2— r2n+l Sm((/3 + ip2n+l)- 



sm K 



(84) 



Let us concentrate on this last expression, for simplicity. 



Since Izr, 



2 

sm K 



sin riK, then if we write Q{n) 



Qi{n){QL^QR)^2e2{n)./Qj:e^, thenO< < 1, 

which implies, in particular, positivity of g. For small k 
it turns out that the second term is negligible (Figure |^). 
We can easily estimate the behavior of the probability 
density profile in this case. For large L and small k we 
have 



Qk{n) 



1 



cos kn/L kn{2n 
■ cos ■ 



1^ 



k^TT^/L^ 



L 



1 + cos ' ^' 1 (pL + qr), 

(85) 

which for fc = 1 can also be approximated as g{n) « 
(1 + ^ sin f (n + i)) (gi + (?h). Figure | shows the ap- 
proximate solution (crosses) and the full solution (di- 
amonds) as well as the probability density of the bot- 
tom states (boxes) and the top states (circles). For the 
smallest resonance (fc = 1) the probability distribution 
achieves maximum around n ^ L/2 where it is approxi- 
mately 2L'^/n^. 

These results are clearly connected to the slow proba- 
bility escape oc To understand them consider a 
plane wave coming from the left with a resonant fre- 
quency going through the open quantum multibaker. 
Thus at every time step we inject the same density inside. 
The wave travels ballistically inside and when it reaches 
the end is mostly reflected, partially transmitted. Due to 
the fast motion inside and slow decay the density accu- 
mulates in the multibaker and reaches the steady state 
when the escape on the right balanced the injection on 
the left. The probability density of the resulting standing 
wave is given by Eq. (p^). 

This result is very striking in comparison with the clas- 
sical case: in the classical multibaker one obtains Fick's 
behavior |, |, § — there is a linear profile of prob- 
ability density. This is also what happens for partially 
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FIG. 8: Profile of the probability density in the steady state 
for the smallest and the fourth resonances k = kn/L with k = 
1,4. Open boxes show the approximate solution, stars stand 
for the full solution. Also shown is the probability density of 
the bottom (full boxes) and the top states (open circles). We 
took Ql ~ 0.1, Qr = 0.9, 77 = 0. The horizontal axis range is 
[-5,105]. 



integrated classical dynamics which we considered in see- 
In particular, the probability density at any 



lie 



tion 

point inside the multibaker is between the densities of 
the reservoirs. 

We defer the complete discussion of the steady state 
solutions to a further work where it will be considered 
together with the semi-classical case in the context of 
transport Here let us only mention that the ap- 

proach to the steady state can be conveniently studied 
as a spectral problem: The evolution equations for the 
quantum multibaker with two waves scattering from left 
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and right can be written as 



(86) 



where := e^''^*^(t), is the matrix representa- 
tion of the open muhibaker propagator following from 
the equations 



foe~ 









Joe 



%in~l,t) 
+ l,t) 
(87) 

and $0 denotes the steady state boundary conditions: 
$0 = [*6(0), $t(0), . . . , ML-1), $t(L-l)]^, $6(0) = A, 
<^t{L - 1) = De''''(^~i\ $b,t(n) = otherwise. The 
solution to this simple afhne problem is 



where Afc are the w-dependent eigenvalues of M^, and the 
\ipk) are the corresponding eigenvectors. In particular, if 
at time the system is empty ^'(0) = 0, then the solution 
is 



(89) 



The steady state is the time invariant part of the above 
solution 



(90) 



The approach to the steady state is given by the eigenval- 
ues of the open multibaker (49), thus it is as slow as the 
escape of probability density, which is consistent with the 
accumulation of large probability density in the system. 
Note that the jiistribution of the absolute values of the 
eigenvalues of is u independent, yet the steady state 
solution does depend on uj. 



THE CLOSED RANDOM QUANTUM 
MULTIBAKER FOR TV = 2 



In this section we extend our discussion of the quantum 
multibaker map for iV = 2 by considering the case where 
the phases, ^q.p, defining the map vary randomly from 
cell to cell. As expected, the random case differs consid- 
erably from the uniform case, since the randomness of the 
phases acts as a disordering mechanism producing a lo- 
calization of the wave function. Unlike the uniform case. 



J 



there is little that can be done analytically for the ran- 
dom case, other than making use of some known results 
for the properties of products of random 2x2 matrices 
^l| , which in this case are only of limited utility. For this 
reason we limit ourselves to a presentation of the results 
of numerical studies. 

The random quantum multibaker map is, for the case 
iV = 2, is defined by the equations 

^fi{n,t+l) = goo(n)/o(n-l)*i(n-l,i) 

+.9oi(ri)/o(n + l)*r(n + l,i), 
*,(n,i+l) = gioin)fo{n-l)^iin-l,t) 

+9ii{n)fo{n + l)-^rin + l,t), 

where the phases in each of the cells are drawn randomly 
from some distribution. Here we use a uniform distribu- 
tion of phases in the unit interval. 

The numerically obtained quasi-energy spectrum is il- 
lustrated in Figure ^ and can be compared with that for 
the uniform case. The quasi-energies associated with the 
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FIG. 9: Spectrum of a the quantum multibaker map for L = 
101 cells. 

eigenstates *(i) = [^'/(O), ^-^(0), . . . , */(L - 1), *^(L - 
1)], are determined by the solution of the following eigen- 
value equation. 







_ 4'r(n) 





5oo("') goiin) 
9io{n) gii{n) 



/o(«- 
fo{n - 



1)^-^(71+1) 



which determines the eigenvalue A. It is interesting to 
note that this equation can be put into a form which 
is reminiscent of the Anderson model for localization 
(see appendix ^). If we define "^{k) by 5';(fc) := 
/o(fc)\E';(fc), 4'r(fc) := /o(fc)^r(fc), we can obtain a set 
of equations that define a generalized Anderson model: 



Join) 



/o(n + l)goi(n) 
A 



$/(n) = <?oo(n)(?i*o(n + l)^i{n - 1) - <?oo(" + l)goi{n)di{n + 1). (92) 
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A similar equation holds for ^'^(ri + 1). We rewrite this 
equation so that it takes the form a dynamical problem, 
where the cell index n plays the role of the time step. 
That is, $i(n+l) = 2\/2e*'^i sm{ip3)%i{n)-e"^-'^i{n-l), 
or, using transfer matrices, 



2\/2e*'^i sin((^3) -e'"^^ 
1 



where the transfer matrix can be written as 



e"f^ 
1 



2\/2sin((^3) -1 
1 



1 

e'^P^ 



$,(n-l 
(93) 



(94) 



The phases are given by 



(pi = {Tr/2)[-~l + (pg{n)ipp{n) - ipg{n + l)ipp{n + l) 

-(Pq{n) ~ (pp{n+ 1)], 
(P2 = -Tr{(pq{n) +ipp{n + l) + 1), 
(fis = K+ {Tr/2)[ipq{n)(pp{n) +ipq{n + l)ipp{n + l) 

-ipq{n) - (pp{n+l)]. 

It can be seen that the eigenstates are localized but 
the localization does not seem to be purely exponential 
for finite L, as illustrated in Figure [l0|, where some states 
are localized over some tens of cells, while others are lo- 
calized over several times as many cells. In longer chains 
the exponential decay of the eigenfunctions is more pro- 
nounced. The exponential decay of the wavefunction far 



A simple estimate can be obtained as fol- 
lows: From Eq. (||) we get |$i(n -I- 1)| = 
|2^/2e*('^l-'^^) sin((^3)$;(n)-|-$/(n-l)|. Thus, on the av- 
erage we have (|$i(n+l)P) = 4(|$;(n)p) + (|$;(n-l)p). 
Therefore, starting from almost every initial conditions, 
on the average we should observe growth of |5'i(n)p 

given by |$i(n + l)|/|$i(n)| « V^ + V^ ~ 2.06. There- 
fore the inverse localization length is approximately 0.72. 
Of course, rather than calculate the logarithm of the 
average we should calculate the average of the logarithm 
but the obtained value is not far off the numerically 
obtained average, which is 0.51 for a chain of 301 cells, 
and 0.56 for a chain of 1201 cells. Figure |Tl| shows the 
distribution of the numbers | In | | over the range 
of the wavefunction where it was appreciably different 
from 0, over all of the eigenstates for a given realization 
of disorder. 
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FIG. 10; Four examples out of 202 eigenstates of a realization 
of the quantum random multibaker of length L = 101 with pe- 
riodic boundary conditions. Absolute values of Vl/;(n), 'l'r.(n) 
are shown. Notice that |^';(n)| — |^'r(?^ 4- 1) | which results in 
particular in the twin peak structures discussed in text. 

from its peak is usually characterized by the inverse lo- 
calization length 



lim In ■ 

|n|— >oo 



+ 11 



(95) 



FIG. 11: Distribution of the ratios A ■- In for a 

I l*i{n)l I 

set of 2402 eigenstates of a realisation of a quantum random 
multibaker of length 1201. 



This distribution reflects the generally broad distribu- 
tions associated with the properties of localized states. 
The difference between the estimated value of the rate 
of growth with the average obtained from the numeri- 
cal distribution is due to the contributions from regions 
where the variation in amplitude from cell to cell is not 
exponential (compare with Fig. [To| ). 

Next, we mention an interesting phenomenon. Wc first 
note that the equations connecting \E'i(n) with ^'/(ri — 1) 
and '^i{n + 1) involve the same phases and are of the 
same form as the equation connecting ^'^(n -I- 1) with 
'i'r{n) and 'i>r{n + 2). When the transfer matrices are 
given in the form of Eq. ( p^ ) one can show that the two 
cases differ only in two random phases, ipi,ip2- However, 
since the equations for ^;(ri) and 'i'r{n + 1) separate and 
have different boundary conditions, there is a priori no 
connection implied between the solutions to the above 
sets of equations. Indeed, after solving the equations nu- 
merically for small L we find that there is no connection 
between them. Thus it may come as a surprise that for 
large L, |\l/;(n)| = |\l/r(?T- -f- 1)|. This is a numerical result 
which accounts for the double peaks in Figure M. To il- 
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lustrate this equality, we show several example fragments 
of eigenstates in Fig. The equality illustrated here is 
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270 280 290 300 200 210 220 230 240 250 



n n 

FIG. 12: The figures show plots of both \<S>i(n)\ and |*r(n+l)| 
superimposed. Fragments of several eigenstates of various 
realisations of quantum random multibakers were chosen {L — 
101 in the first figure, L = 301 in the following figures; the 
parts not shown are zero within the numerical precision). 

satisfied to within all the allowed precision for L > 100, 
and is independent of the shape of the eigenfunction. 

The physical explanation is actually quite simple ^T\ . 
The probability current through the boundary between 
cells n and n + 1 is given by J„|„+i = — 

|4'r(«+l)p. Therefore the difference between |^/(n)| and 
|5'r(?^ + 1)1 implies non- vanishing current. While such 
currents can arise in principle one expects they should 
be negligible in disordered systems, which leads to the 
equality which we observe numerically. 



VI. SUMMARY OF RESULTS AND 
DISCUSSION 

With this paper we have begun our study of transport 
properties for quantum multi-baker maps. Our central 
motivation for this study is, of course, to use a sim- 
ple, classically chaotic system with transport and en- 
tropy production and to explore the differences between 
the classical system, and its quantum counterpart. The 
multi-baker map provides a convenient model for this 
study. 

In order to set the stage for further work, we consid- 
ered here the extreme quantum limit of the QMB where 
each unit cell contains just two quantum states, in ei- 
ther the position or momentum representation. Further 
work will now be devoted to the properties of the QMB 
for smaller values of h, including a study of the semi- 
classical regime where the Planck constant is very small, 
and the transition to the classical limit. 

In order to explore a variety of versions of the QMB, we 



used the fact that the Weyl quantization procedure for 
maps on torus allows some freedom in the choice of phases 
for the quantum states. This freedom can be thought of 
as a freedom in the location of the quantum states in the 
position and momentum representations, or as the effects 
of Bohm-Aharonov current loops on the wave functions 
in each unit cell. For large h and for the uniform QMB, 
we have been able to proceed with analytical calculations 
for the quasi-energies, and for the quantum states for dif- 
ferent boundary conditions. These calculations showed 
ballistic transport of probability amplitudes within the 
chain, due to the Bloch wave structure of the eigenstates, 
as well as a slow decay of probability from an open, finite 
chain, due to the difficulty that long wavelength modes 
have in escaping through the exit channels. The random 
phase case exhibits localization of wave functions and 
the appropriate equations were shown to be close to the 
basic equations of the Anderson model of localization. 
All of these results show that the behavior of the QMB 
for = 2 is quite different from that of the classical 
multi-baker map, which exhibits normal, diffusive trans- 
port. For larger N, smaller h, we should see behavior 
that more closely approximates the classical behavior, at 
least on the logarithmic time scale. 

There are quite a large number of directions for future 
studies of the QMB. For the case of ft, = 1/2 we have stud- 
ied only the uniform and random phase cases. One can 
also consider models with periodic distribution of phases 
over several cells, models in which there is a well defined, 
systematic progression of phases from one cell to the next, 
or where the phases in the cells are incommensurate (os- 
cillating with irrational periods) . The situation for larger 
A^ is such that there are A^ quantum states in each cell, so 
that there will then be several transport channels in the 
QMB. It remains to be seen what the properties of these 
systems will be, both as functions of N as well as func- 
tions of the phases of the wave functions in the cells. As 
h approaches zero, one can study the semi-classical limit 
and the approach to classical, chaotic behavior of the 
QMB. Clearly, we should obtain the diffusive behavior 
independently of the distribution of phases, that is both 
for the uniform (ballistic) case, as well as for the random 
(localized) case. Such a study should reveal for this sim- 
ple system whether the diffusive behavior arises directly 
through the semi-classical limit, or whether one needs to 
assume additional mechanism — such as interaction with 
environment, or decoherence — to regain classical prop- 
erties. The QMB is simple enough to consider the clasical 
limit for uniform phase, random phase, and incommen- 
surate phase model, among others. The investigation of 
these questions will be the subject of further papers in 
this series. 
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APPENDIX A: EIGENVALUES OF THE OPEN 
QUANTUM MULTIBAKER LIE INSIDE THE 
UNIT CIRCLE 

Due to the escape of probability density, the eigenval- 
ues that determine the time dependence of the probabil- 
ity density in each cell move to the interior of the unit 
circle. A simple proof of this fact can be given as follows. 
Let = 0) be a normalized eigenstate of the open 
quantum multibaker. Since the multibaker dynamics re- 
quires that *f,(n-l,i-|-l) = fo^i{n,t), $t(n-|-l,< + l) = 
fa'i'rin, t), for n = 1, . . . , L — 2, it follows that 

|*b(n-l,i+l)|2+|*,(n+l,t+l)|2 = \^i{n,t)\^+\^,{n,t)\' 

(Al) 

The probability of the system being in cell n is given 
by gin,t) \'ifi{n,t)\^ + \^r{n,t)\^ = \^b{n,t)\^ + 
|5't(n, where we have used the unitarity of the trans- 
formation G2^- In the boundary cells we have 

|*f,(0,i+l)|2 + |$t(i^-l,f+l)|2 = 0, 

|M',(l,i+l)|2 = \^l{0,t)\^ 

\ML-2,t+l)\^ = |vI/,(L-l,t)|2. 

(A2) 



By adding up the equations (Al) for cells n 



1, 2, . . . , L - 2 and (|A3) we obtain 



L-1 



L-1 



J2 0{n,t+l) = J2 Q{n,t)-{\^r{0,t)\^ + \'^i{L-l,t)\^). 

(A3) 



n=0 n=0 

But is a normalized eigenstate at t = 0, therefore 



L-1 



Q{n,l) = |Ap = 1 - (1*^(0, 0)p + l^iiL - 1,0)1^ 



This implies 1- |A|2 = (|*^(0, 0)|2 -h l^-^i - 1, 0)P) > 0. 
Thus < |Ap < 1. Suppose now |A|2 = 1. Then 
*r(0, 0) = *;(i- 1, 0) = 0. Thus, from Eq. (|l|) and (|2|) 
= *6(0,0) = .9oi*i(0,0) ^ *K0,0) = 

^^(0,0) = 0. It follows that also *t(0,0) = 0. But 
$t(0, 0) = A"i$t(0, 1) = A-i/o*r(l, 0). Thus ^',(1, 0) = 
0. Also, *6(1,1) = /o*i(0,0) = = A.goi*i(l,0). There- 
fore ^'/(1,0) = *r(l,0) = 0, and so on. Thus the as- 
sumption |Ap = 1 leads to the eigenstate being iden- 
tically which cannot be normalized. Therefore all of 
the eigenvalues lie inside the unit circle < |Ap < 1. 
It is easy to identify the kernel of Ml since A = im- 
plies |*^(0,0)p |4',(L - l,0)p = 1. Thus the kernel is 
spanned by vectors with everywhere apart from 5'r(0) 
and *;(£ - 1). 

We are thus led to the conclusion that of the 2L eigen- 
states exactly two span the kernel, and the eigenvalues 
corresponding to the remaining 2L— 2 satisfy < |A| < 1. 
Note that the above arguments are independent of the 
distribution of phases and thus apply also to the open 
random quantum multibaker. Generalization to arbi- 
trary is obvious. 



APPENDIX B: ANDERSON MODEL 



Consider a one-dimensional Schrodinger equation on a 
lattice of lattice constant a lP2i : 



2ma' 



■[^'„+i + *„_i-2«'„]-|-t/„*„. (Bl) 



A time independent equation can be written as {E — 
T4)*„ = ^'„+i-|-*„_i, where E = 2 ^ 2ma^E /h^ , V„ = 
2ma^Vn / fi^ ■ We can also write it using transfer matrices 



(B2) 



which has very similar form to the equation (P3[). A 
general discussion of one-dimensional disordered models 
can be found in 131], p3l 
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